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This report outlines the theory used in FAS_TER-III, a. Monte Carlo 
computer program for the transport of neutrons and gamma rays in com- 
plex geometries. The code has the additional capability of calculating 
the minimum weight layered unit shield configuration which will meet a 
specified dose rate constraint. It includes the treatment of geometric 
regions bounded by quadratic and quadric surfaces with multiple radia- 
tion sources which have a specified space, angle, and energy dependence. 
The program calculates, using importance sampling, the resulting number 
and energy fluxes at specified point, surface, and volume detectors. 

Results are presented for sample problems involving primary neu- 
tron and both primary and secondary photon transport in a spherical 
reactor -shield configuration. These results include the optimization of 
the shield configuration. 


Section 1 

IHTRODaCTION AND SUMMARY 

The original FASTER program (ref. l) contained 
a number of new techniques which provided the capa- 
bility of obtaining accurate radiation levels at 
specified points in complex geometries. Prior use 
of FASTER indicated a need to broaden the overall 
program capabilities, automate the importance 
sampling, increase the computational efficiency, 
and revise the users manual. This revised program 
has been designated FASTER-III to distinguish it 
from earlier versions. 

A specific program capability permitting the 
calculation of minimum weight layered unit shield 
configurations for mobile nuclear reactor applica- 
tions, e.g., nuclear propulsion for aircraft, sur- 
face effect vehicles, and spacecraft has recently 
been developed. The basic Monte Carlo tranf|)ort 
method was extended to Include a calculatloh of 
partial derivatives of the radiation fluxes with , 
respect to specified shield dimensions. These de- 
rivatives are then used to define exponential re- 
lationships used in the shield optimization pro- 
cedure. This optional program feature is described 
more completely in Section 2. 

Data preparation is simple, with very little 
judgment required to set up the importance sampling j 
for most problems. The code also has a \mit shield i 
weight optimization capability. 

Particularly noteworthy features of FASTER- 
III are the following: 

(1) A calculation of optimal importance sam- 
pling parameters based on partial derivatives of 
the variance (Section 2.3). 

(2) The acceptance of data in either fixed or 
variable field formats including the ANISH-DTF for- 
mat for neutron cross sections. 

(3) The calculation of time -dependent neutron ; 
and photon transport (using time moments and/or 
time intervals) including an optional exponential 
atmosphere. 

(4) The improvement and addition of importance 
sampling models with the various importance sam- 
pling parameters built into the program. 

Various program features are described in 
Refs. 2 to 6. 

The application of the FASTER-III program to a 


shield optimization problem is discussed in Section 
3. The problem involved a spherical reactor-shield 
configuration and included primary neutrons and 
both primary and secondary photons. Conclusions 
and recommendations are presented in Section 4. 

Section 2 

ANALYSIS 

The techniques used in calculating optimum 
shield configurations and optimxim importance sam- 
pling parameters are summarized helow. The dis- 
cussion is given in three parts: dose rate deriv- 

atives with respect to shield layer thlc3snesses, 
optimization procedures, importance parameter op- 
timization. 

2.1 Dose Rate Derivatives 

The dose rate at a point detector y for a 
specified reactor shield configuration is written 
as: 




( 1 ) 


j = l 


where J is the total number of energy groups for 
both neutrons and photons (including secondaries), 
Tj(z) is "the particle flux in the energy group, 
and Rj is the response function to convert from 
flux to dose rate. The rate of change of the dose 
rate with respect to a shield layer thielsness is 
simply 
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where L is the total number of shield layers and 
tx is the thickness of the Xtt layer. The equa- 
tion used by the program for determining the flux 
is written as : 
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where N is the total nimiber of histories tracked 
via the Monte Carlo method, k is the number of 
particle collisions, is the position of the 

kth collision of the n th history, sJjoj(ujqi) the 
number of particles in the energy group 

emerging from zj^ in the direction of the 

detector per unit solid angle, and rep- 

resents the material and geometric attenuation 
kernel for particles in the energy group 

going from to the detector. 

The partial derivative of the flux with re- 
spect to the shield layer thickness is 

simply: 
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The summations are a minor part of the calculation. 
Therefore, the notation is simplified by concentra- 
ting on the elements in the summation 

^ [3!kn("kn)Kj(am^ z)] (5) 


where Sjkn represents the contribution to the 
flux in the energy group from the k th col- 

lision of the n th history. This equation is re- 
written as 
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The second term in brackets involves the attenu- 
ation kernel 
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where M is the total number of regions traversed 
from Zj^ to the detector, s is the path length 
for the mth region traversed, is the total 

cross section of this region for particles in the 
J th energy group, and s is the total distance 
from zjjjj^ to the detector, i.e.. 
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A substitution of this kernel gives: 
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The partial derivative of the partial path 
length S]u with respect to the shield layer thick- 
ness tj is zero unless the mth region tra- ' 
versed is affected by a change in t^. In partic- 
ular, if t^ is a characteristic dimension of the 
region, i.e., its thickness, then 
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where is the cosine of the angle measured 

from the surface normal nbmp with which the par- 
ticle crosses the boundary of the region. 

In the strict sense, the change of the thick- 
ness of one shield region can affect other shield 
regions. In particular, for a spherically symmet- 
ric reactor -shield configuration, an increase in 
the thickness of a shield region forces a movement 
of all shield regions having a larger radius. The 
inclusion of these effects in the above equation 
unnecessarily complicates the analysis and the cal- 
culations. The primary effect of changing a shield 
region dimension is to change the number of mean 
free paths which particles have to traverse in 
reaching the detector. Therefore, in calculating 
the derivatives, only the effect of the material 
attenuation is treated. 

The derivatives at a specific boundary cross- 
ing m' then simplify to: 
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where m' is the index of a region having tj as 
a dimension. The partial derivatives of the par- 
ticle weight with respect to the shield dimensions - 
the first term in brackets in Eq. (6) - are zero 
at the point of origin of all primary particles. 

For subsequent particle collisions, the deriva- 
tives are calculated using the relationship be- 
tween particle weights on subsequent collisions: 
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where S* jj_3_ ^(Zkn) number of particles 

coming ou't of ^ the previous collision point in the 
direction Vj^j^ and in the i^ energy group, 
^i(2k-l ni 5kn) attenuation kernel between 

particle collision points, Tij(zjmj 3m ' Hkn) is 
the scattering kernel for transfer of particles from 
group i to group j, and l^(zioi) is ills probability 
density function used in selecting the collision 
point. 


A straightforward substitution g^es 
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"After some manipulation, this reduces to 
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where 


^IJkn 
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The first term in brackets in Eq,. (14) is the 
same partial derivative for collision k-1 as the 
partial derivative now being calculated for colli- 
sion k. Therefore, it is known, either identi- 
cally zero for k=0, or as determined from Eq. (14) 
for k > 0. The second term in brackets in Eq. (14) 
is similar to the second term in brackets in Eq. (6) 
and is therefore determined by Eq. (ll). The last 
term in brackets involves the definition of the 
probability density function used to select the 
collision point Zj^. 

The probability density function for a colli- 
sion point has the form 
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where <l5n(Sm) ^ probability density function 
used to select the particle direction, s = |zj^ - 
2k -1 nl distance of the selected collision 

poin"!; from the previous collision point, A(s) Is an 
importance factor for each region which changes 
dlseontinuously at region boundaries, and a(s) is 
an effective cross section which changes discontln- 
uously at region boundaries and which may change 
continuously within a region. 

The derivative of the logarithm of P^(z^) 
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involves only those terms which change when a 
shield dimension changes, i.e.. 



Let sj denote the distance to a boundary in- 
volving the shield dimension. If the first 

term on the left side of Eq. (l?) is affected by a 
change in this shield dimension, i.e. if s > sj, 
then 

T- r®a(s')ds' = -a(sj) ^ 

Stj L *X) J 


= -a(s^) (18) 
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where a(s^) is the effective cross section at the 
boundary of the shield and (Jjjsn is the cosine the 
particle path makes with the outer shield normal. 

If there is any crossing involving the Ztti shield 
dimension, the second term in Eq. (l8) will always 
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Curved shield surfaces may he crossed more than once 
along the path between two particle collision 
points. Therefore, a summation of Eqs. (l8) and 
(19) over every intersection involving the I th 
shield dimension is required to completely evaluate 
Eq, (17). 

2.2 Optimization Procedures 



The shield optimization calculation yields the 
set of shield layer thicknesses t' = (t4, t^, 
t), t£) such that the dose rate, D(t'*‘), meets 

the dose constraint. The Monte Carlo calculation is 
performed for an initial set of shield layer thick- 
nesses t = (tj, t 2 , t^, ..., t]^) and yields a 

set of fluxes, cpj ( t) , J = 1, ?, . . . , J and deriva- 
tives, Scpj(t)/St^, j = 1, 2, ..., J; 7 = 1, 2, ..., 
L. The assumption is made that the fluxes vary 
exponentially with respect to shield dimension 
changes in the form 



0=1 
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(26) 


cpj(t') = cpj(t) exp^j • (t' - t^ (20) 

where aj = (ajj^, a^g, ..., Sjj^). It follows that 

9j(t) exp[aj . (t- - t)] ^[aj . (t ' -t)] 


and for spherically symmetric shield: 



(27) 


= cpj(t')ajj (21) 

In particular 



The weight is also expressed as a function of 
the shield layer thicknesses. The weight is de- 
noted hy tf(t') and for spherically symmetric 
shields : 



where is the density of the 7th shield re- 

gion and rp is the, minimum shield radius. 

The purpose of the optimization procedure is 
to minimize the weight W(t ' ) subject to the dose 
rate constraint D(t') = Dp where Dp is a speci- 
fied dose rate. At this optimum, a small weight 
perturbation in any layer causes the same dose rate 
change. The rate at which dose rate changes with 
respect to a shield weight change in the Jth layer 
is given by 


In arriving at the optimum shield, the total 
shield weight is built up in increments of weight 
in. Each increment in shield weight is always 
associated with a particular shield layer thickness. 
At each iteration, the particular shield dimension 
is selected by examining the values of the shield 
weight quality factors, Q,- Each factor Q^, rep- 
resents the approximate change in dose rate per 
unit change in weight corresponding to a change in 
the 7to shield dimension. Negative Qj ' s are the 
most usual and correspond to shields for which an 
increase in weight - and shield dimensions - gives 
a decrease in dose rate. Positive Q^'s can occur, 
however, and correspond to shields for which an in- 
crease in weight also increases the dose rate. 

If, at a particular iteration, the dose rate 
is above the dose rate constraint, the minimum 
shield weight increment would correspond to the 
least positive value of those Qj ' s for which 
0,1 > 0 and for which tj > t^(min), where tj(min) 
is the minimum value of the 7th shield layer 
thickness. If such a Q, exists, the dose rate 
can be decreased while also decreasing the shield 
weight the maximum amount. If there isn't such a 
Oj, the next best procedure is to find the most 
negative of the Oi ' s for which Qj < 0 and for 
which tj < tj(maxJ, where t^(max) is the maximimi 
value of the 7;tt shield layer thickness. A 
change in that Qj would give the maximum decrease 
in dose rate per unit increase in weight. 

If the dose rate is below the specified dose 
rate at a particular iteration, the minimum shield 
weight increment would correspond to the least 
negative of those Qj ' s for which < 0 and for 
which t| > t^(min). If such a Qj exists, the 
dose rate can be increased while decreasing the 
shield weight the maximum amount. If there isn't 
such a Qj, the next best procedure is to find the 
most positive of those ' s for which 0i > 0 
and for which tj < t j (max) . A change in that 
would give the maximum increase in dose rate per 
unit increase in weight. 


4 



Assuming a particular value of the Q-j's 

is selected tlirough the above arguments, the cor- ' 

responding shield dimension t^J^ is changed by a 

maximum amount At where At„ is calculated as 
m m 


At 




aw(t ' ) 

St' 


(28) 


If this change would put outside one of 

its specified limits, the value of t^J^ would be 
set to that limit, i.e., tj„(min) ^ ^ tj^(max). 

The shield weight increment AW is calculated as 
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subject to the constraint that |ZW| < where 
ZWq is a specified maximum shield weight incre- 
ment per iteration. Note that AW, and therefore 
Atjj^, may be positive or negative depending on the 
value of and whether the dose rate is above 

or below the dose rate constraint. 

Once a shield layer thickness is changed, the 
dose,- weight, and their derivatives are re- 
evaluated and the entire process is repeated. The 
optimization would he discontinued in several 
ways. If the dose rate equals the dose rate con- 
straint within the relative error of the original 
Monte Carlo dose rate calculation, the program will 
proceed to the next problem - which may he identi- 
cal except with more histories to tighten the con- 
vergence of Monte Carlo calculations. Similarly, 
if all shield layer thicknesses have reached their 
minimum or maximum values, and if the optimum 
shield cannot he determined with these constraints, 
the program would again proceed to the next prob- 
lem. Finally, if the dose rate and dose rate con- 
straint are decades apart in value, the program 
would reevaluate the fluxes and their derivatives 
hy Monte Carlo every time the dose rate changed by 
more than a specified factor during the optimiza- 
tion procedure. 
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Thus the partial derivative of the relative 
error with respect to the parameter a is: 
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The dose rate from the n th history is given hy 
J 


j=l k 


(34) 


where J is the total number of energy groups, k 
is the number of particle collisions, Rj is the 
flux to dose rate conversion factor for^the j th 
energy group, and is the flux in the j th 

group from the k th ‘^collision of the n th his- 
tory. Since 
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2.3 Importance Parameter Optimization 

The optimization of the Importance sampling 
must he performed for some function, e.g., dose 
rate, of the energy-dependent fluxes since there 
is a different optimum for every Initial particle 
energy. Therefore, assume that a minimum variance 
calculation of the dose rate is required where 


the calculations required to evaluate Eq. (33) all 
involve the summation of terns which involve 



% = iJ'Dn (30) 

n=l 

where W is the total number of histories and Djj 
is the dose rate from the n th history and Dj} 
is the average value of the dose rate after K 
histories. The relative error of this dose rate is 
given hy 


The remainder of the analysis, therefore, can 
he concentrated on the partial derivatives of the 
fluxes. All other operations which must he per- 
formed are given above. 

The fluxes typically depend on the detector 
position 2 ^ so the equation for the particle flux 
is written as 

'Pjkn(z) = sW22kn)Kj(Sfcn^ £> (37) 



Taking the logarithm of this equation and 
then performing a formal calculation of the partial 
derivative with respect to an unspecified param- 
eter a yields 


The transport kernel Kj(zjj^, 2 ) does not in- 
volve any importance sampling parameters so that 


^jkn^^^ _ ^^jhn(Hkn^^j (.Skn^ Z) 
3a. Ba 
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This equation can also he written as 

^ in S]j^(ui^) (39) 
oa oa 

Without going into great detailj it turns out 
that the particle .weight Sjj^j(uj^) is composed of 
a purely analytical numerator, ^lkn(likn) 
nominator which is the product or all the prohahil- 
Ity density functions used to select the collision 
points, i.e.. 


a* /„ ^ _ ^,1kn^-kn^ 

^Jkn'-Hkn^ 
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Therefore, 
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n ^tn^%n^ 
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Since Vj]j^(ujjjj) does not explicitly Involve 
any Importance parameters, it follows that 


to a set of parameters which can have a reasonably 
simple role. These parameters consist of the rel- 
ative importance Ij, of each region. Normally 
these parameters are all equal. However, in asym- 
metric problems, it turns out that some regions are 
much more important in terms of their scattering 
contributions to a detector. Therefore, these im- 
portant regions have a larger value of I . 

The region importance enters into tSe selec- 
tion of a collision point through the following 
probability density function: 
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where r is the region in which the collision 
occurs (selected at random), p*(s) is the piecewise 
continuous probability density function in this ’re- 
gion at the selected collision point (a distance s 
from the previous collision point), H is the total 
number of regions in which the collision could have 
occurred, and is the integral of p^( s ' ) over 

the partial path length in region h. 

Calculating the logarithm of each side of the 
equation yields: 
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■ The partial derivative of Eq. (47) with re- 
spect to the specific importance parameter Ig - 
the relative importance of region g - yields 


Therefore, Eq. (39) can be re-written as 
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Moreover, the partial derivatives are energy- 
independent so that Eq. (36) becomes 
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The evaluation of the partial derivatives of 
the probability density functions can be written 
as 
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where 6gjj =0 if region h is not region g and 
5gg = 1- 

Thus Eq. (48) is evaluated during the random 
selection of the k th collision point and the 
final term necessary to evaluate Eq. (45) and all 
preceding equations has been determined. ■ 

The above analysis is used to calculate the 
partial derivatives of the relative error of the 
dose rate with respect to the relative importance 
Ij. of each geometric region, and a similar analy- 
sis is performed for the other importance sampling 
parameters. The result of the complete Monte Carlo 
calculation is a set of partial derivatives which, 
for the region importance, are given by 
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Sa 
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At the k th collision, the first term on the 
left side of Eq. (45) is known, identically zero 
if k = 0. Therefore, the analysis is completed 
after examining the calculation of the second 
term. 

At this point it is necessary to Identify the 
particular importance parameter a. Since most of 
the importance sampling parameters have fairly 
Involved roles, the technique will be applied here 
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(49) 

where SDj^/SIj, Is obtained from Eq. (44) using 
Eqs. ( 45 ) and (46). 

After the calculation is completed, optimal 
values of the importance sampling parameters are 
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calculated ty requiring that the relative error he 
zero - not actually achieved of course. 

By a first order expansion 



where R is the total number of regions. A simple 
gradient analysis says that should he 

proportional to SEjj/SI^ so that “ 


I, + C 
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(51) 


where, hy substitution into Eq. (50), 


C = 



where z is a general phase 
space coordinate 


(5K) 


The program prints the optimum values of 
and other importance parameters after completing— 
the Monte Carlo flux calculation. This analysis 
is performed for every response function. After 
more experience is obtained with the technique, the 
program could be modified to change these param- 
eters internally corresponding to a specified re- 
sponse function. 


Section 3 

SAMPLE ERCBLEM RESULTS 


Two problems were Investigated using the 
shield optimization capabilities of the FASTER-III 
program. Both problems involved a spherical 
reactor-shield configuration and included primary 
neutrons and both primary and secondary photons. 

The two problems were similar except for the 
power level, 375 MW and 600 MW respectively. Both 
problems used a flat radial distribution for the 
primary neutron and photon source distribution. 

The primary photon source included an infinite 
operation equilibrium fission product term. 

The core radii for the two problems were 82.38 
and 96.38 cm respectively, corresponding to a power 
density of 4.53 MW/ft^. Following the core was a 
7.62 cm Be reflector; a 5 cm depleted uranium 
shield; three depleted uranium-borated water shield 
layers of 57, 15, and 15 cm thickness and 6.4, 

4.6, and 2.8 gm/cm^ density respectively; and a 
117 cm borated water shield. This base line shield 
configuration was based on parameters obtained from 
SANE-SAGE calculations and subsequent calculations 
using the UNAMIT program. Ref. 7. The reactor- 
shield compositions are given in Table 1. 

The primary neutron transport calculation 
utilized multigroup cross sections for 26 energy 
groups. Fifteen energy groups were utilized for 
both primary and secondary photons. The secondary 
production cross sections included both inelastic 
and capture gammas. 

These initial configurations were each analy- 
zed for a point detector 30 ft from the core cen- 
ter by following approximately 500 energy- 
dependent packets of primary neutrons and photons 
and approximately 7000 packets of secondary pho- 
tons. The dose rates obtained from these calcu- 


lations are tabulated in Table 2 including a break- 
down by secondary source region. Each of these 
problems required about 28 minutes on the UNIVAC 
1108 computer. 

The basic calculated dose rates and dose rate 
derivatives were also used by the FASTER-III pro- 
gram' to calculate the minimum weight shield con- 
figuration which would give a dose rate of 0.25 mr/ 
hr at the specified detector point. The final ■ 
shield configurations following the optimization 
are given in Table 3. 

In both cases, the optimum shield configura- 
tion is significantly different from the base line 
configuration. Since the base line configuration 
was not generated by the FASTER-III program it is 
difficult to discuss many factors entering into 
that calculation which would account for the dif- 
ferent optimal configuration. It is noted, how- 
ever, that the . base line configuration was gener- 
ated using parameters corresponding to a calculated 
dose rate an order of magnitude below the specified 
dose rate constraint. Ref. 8. As such, the base 
line configuration used in the FASTER-III program 
was determined from an extrapolation of a different 
base line configuration. 

A more critical critique can be made of the 
FASTER-III results independently. First it is 
noted that neither problem saw a significant con- 
tribution, less than a few percent, from photon 
sources in the core region. In fact, the 600 MW 
reactor dose rate from this source was about a fac- 
tor of two less than it was for the 375 MW reactor. 
This difference is ascribed to the problem statis- 
tics since core photon sources see approximately 
30 mean .free paths of shield material. Therefore, 
it is doubtful if this dose rate component is con- 
verged within a factor of two after only 500 pack- 
ets but this does not introduce a significant error 
since the original contribution was only two per- 
cent of the total dose rate. 

The small contribution from core photon 
sources decreases the amount of high Z shields 
required around the core. Therefore, both problems 
gave a significant change in the first two shield 
dimensions during the optimization. In the 375 MW 
problem, the first mixture of depleted uranium- 
borated water (p = 6.4 gm/cm^) was eliminated en- 
tirely. In the 600 MW problem, the depleted ura- 
nium and most of the first mixture were eliminated. 

The main difference between the two FASTER-III 
calculations was the shift in the placement of 
lighter shield mixes towards the core for the 
600 MW problem. An examination of the secondary 
photon dose components indicates that the contri- 
bution from the outer two shields was about 25 per- 
cent for the 375 MW reactor and almost 50 percent 
for the 600 MW reactor. Since these sources de- 
pend on the neutron attenuation through the closer 
regions and since lower effective Z materials 
are better neutron attenuators on a weight basis, 
the 600 MW problem tends to replace high effective 
Z material with a lower effective Z material. 

The differences in the contribution from 
secondary sources in the outer shield regions is 
greater than expected for the nominal difference 
in the core region. Therefore, much of the dif- 
ference in these sources must be ascribed to sta- 
tistical variations. In fact, both problems had 
approximately 25 to 30 percent calculated relative 
error in the total photon dose rate. It should be 
noted that the FASTER-III program includes a num- 
ber of Importance sampling techniques which could 
be used to decrease this error. However, both 
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problems were run using the tuilt-ln definitions of 
importance parameters. Alternatively, more his- 
tories could have heen used although the computer 
time requirements would have hecome excessive. 

Section 4 

CONCLUSIONS AND RECOMMENDATIONS 

The FASTER-Ili program was developed to calcu- 
late neutron and photon fluxes at specified points 
in complex geometries. Alternatively, it can also 
calculate fluxes averaged over specified surfaces 
and volumes. The program was designed such that 
data preparation is simple and so that very little 
judgment is required to set up the importance sam- 
pling for most problems. The FASTER -III program 
satisfies these requirements very well. 

The shield weight optimization capability in- 
cluded in the FASTER-III program permits the cal- 
culation of both base line radiation levels and op- 
timal shield thicknesses all in a single computer 
run. However, the very large attenuation factors 
Involved in the demonstration problems yielded some 
questionable results. In particular, the statisti- 


cal differences in the relative contribution from 
various secondary source regions caused correspond- 
ing variations in the relative distributions of 
shield materials. Of course the statistical varia- 
tions would be less in problems with less overall 
attenuation. 

■ The effect of statistical differences on the 
shield optimization can be reduced by following 
more packets. However, the computer times start to 
get excessive if this is the only approach used. 

It would be more fruitful in terms of the routine 
application of the program to expend some effort 
towards altering the Importance sampling. 

The FASTER-III program has the capability of 
calculating optimal importance parameters based on 
partial derivatives of the variance. This feature 
can be used in determining better Importance sam- 
pling parameters for shield optimization problems. 
In fact, the overall program efficiency could be 
Improved if this feature was utilized on a wide, 
variety of problems with the results being used to 
improve the built-in importance sampling models and 
parameters . 
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8FKBRICAL REftCTOH-SHXEIS CON?IQURATION 



COMPOSmOlB (10^** atems/en^) 


CORE 

OPUCTO 


MIX 1 

MIX 2 

MIX 3 

nmii 




SMISU) 

SHIEID 

■SHTBTJt 

SHIEID 

H 

0.01976 

0.0 

0.0 

0.0451 

0.0516 

0.0580 

0.G645 

B.S 

0.0 

0.120 

0.0 

0.0 

0.0 

0 0 

0.0337 

B 

0.0 

0.0 

0.0 

0.000671 

0.000766 

0.000862 

0.0009S8 

0 

0.0118A 

0.0 

0.0 

0.0226 

0.025S 

0.0290 

O.Q 

Al 

0.0512 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 








18 

Zr 

0.01744 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.000979 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

d“38 

0 . 00007 S 

0.0 

0.0482 

0.01446 

0.00964 

0.00462 

0.0 


i 


TABLE 3 

RESUUD3 OP PA3T15R-III 3HIELD OmMlZATtOK 
(0.25 mr/hr at 3<'' feet) 


RESULTS OF PASTHl-III BASE LINE CALCUIATI<»IS OP REACTCR 
SHtEU) COKFIGURASIOHS AT 30 FEET PRCW .CORE CENTER 


DOSE RATE 
CO1P0HEMT 


DOSE RATE COHTRIBOTIONS AT 30 FEK- 
FR(»I CCSE CEWTER 


375 MW 

REACTOR 

fmr/hr) 


600 MW 

REACTOR 

(mrAir) 


Photon Source Region 


Core 

0.009 

0.004 

Reflector 

3.5x10'® 

6.3x10' 

Depleted Uranium 

3.2x10‘5 

1.3x10* 

Mix 1 Shield 

0.018 

0.026 

Mix 2 Shield 

0.062 

0.075 

Mix 3 Shield 

0.017 

0.063 

Berated Water Shield 

O.Oll 

0.022 



375 MW 
REACTOR 

600 MW 
REACTOR 


Initial 

Pinal 

Initial 

Pinal 

Dose Rate (mr/hr) 





Photon 

0.120 

0.126 

0. 187 

O.L53 

Neutron 

0.020 

0. 124 

0.027 

0.097 

Total 

0.140 

0.250 

0.214 

0.250 

ajleld Weight (lO^kg) 





Depleted U ' 

10.2 

12.6 

13.8 

0.0 

Mix 1 

71.2 

0.0 

89.2 

6,6 

Mix 2 

22.1 

52.4 

26.4 

52.4 

Mix 3 

16.1 

12.2 

19.0 

63.1 

Water 

86.7 

80.3 

2I1I 


Total 

206.3 * 

157.5 

246.1 

207.4 

Shield Thickness (om) 





Depleted U 

5.0 

6.1 

5.0 

0.0 

Mix 1 

57.0 

0.0 

57.0 

7.0 

Mix 2 

15.0 

57.3 

15.0 

48'. 4 

Mix 3 

15.0 1 

13.5 1 

15.0 

51,4 

Water 

LI7.O 

120.8 

117.0 

98.4 


Total Photons 


0.120+0.034 ' 0 . 187 * 0.054 


0.020+0.002 


0.027+0.003 


Total 


0.214 


NASA'Lowis 



